Order a 2 corrections to the decay rate of orthopositronium 
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Order a 2 corrections to the decay rate of orthopositronium are calculated in 
the framework of nonrelativistic QED. The resulting contribution is found to be in 
significant disagreement with one set of experimental measurements, though another 
j experiment is in agreement with theory. 
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The discrepancy between theory and experiment for the decay rate of orthopositronium 
has long been one of the outstanding problems in precision QED. The theory is given by 
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r o _ Ps = r [l + A- + ?-\na + B(-) 2 - ^—\n 2 a] (1) 
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where the lowest order decay rate is given by 

r 2(^-9)^. (2) 

y 7r 

The value of T is 7.211 169 /is" 1 . This is about 2.3 percent above the experimental 
determinations at Ann Arbor of 

r o _ Ps (Gas) = 7.0514(14)/i S - 1 [1] (3) 

r o _p s (Vacuum) = 7.0482(16)/is~ 1 [2] (4) 

and 2.4 percent above the somewhat less precise Tokyo measurement 

r o _ Ps (Si0 2 ) = 7.0398(29)/is" 1 [3]. (5) 

The great bulk of this difference is accounted for by the one- loop correction 0, which has 
been evaluated with high accuracy in M to be 
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A = -10.286 606(10). 



(6) 



Including this -2.39 percent effect along with the logarithmic terms of order a 2 T |J and 
a 3 r 0, which contribute -0.01 percent and -0.0004 percent respectively, gives a decay 
rate of 7.038 202 yus -1 . The remaining difference with the Ann Arbor experiments, with 
the 2.3 percent difference reduced to -0.1 percent, requires a rather large positive value for 
B (339(36) for the gas experiment and 257(41) for the vacuum experiment), which is the 
discrepancy mentioned above. The Tokyo experiment, on the other hand, which disagrees 
with the first two by several standard deviations, is consistent with a small value. While the 
experimental situation clearly requires more work, it is also obviously important to directly 
evaluate the constant B. 

While the need for calculating B has been clear for two decades, there are two difficulties 
that have prevented its evaluation until now. The first is simply the large number, 81, of 
diagrams that contribute, many of which have two-loop ultraviolet infinities. More impor- 
tantly, a number of these diagrams have a serious kind of infrared divergence associated with 
the fact that positronium is a bound state. 

These problems have recently been overcome for the case of parapositronium decay ||. 
The main theoretical tool used in this work is Nonrelativistic Quantum Electrodynamics 
(NRQED) |§. This approach allows the high-energy part of the problem to be treated 
as an on-shell scattering process. In this case the complications of the bound state are 
not present, and Feynman gauge can be used. The low-energy, bound state part can be 
treated in Coulomb gauge with a small set of operators that describe relativistic and QED 
corrections to a Schrodinger problem. A matching procedure carried out with free particle 
scattering amplitudes then allows the two parts to be combined into a complete calculation. 

The present calculation, while similar in spirit to that of the parapositronium calculation, 
regulates long wavelength singularities by giving the photon a small mass m e A: reference || 
instead uses dimensional regularization. We have chosen to use a slight variation of a recent 
NRQED calculation of one-photon annihilation contributions to ground state positronium 



hyperfine splitting |L0| that uses a photon mass, as it is easily generalized to the decay rate 
calculation. 

In NRQED the annihilation of orthopositronium can be accounted for by modifying the 
amplitude for one-photon annihilation, 27ra/m 2 , to 



^,i) = ^(i- 4ia Y" 9) )- m 



The independence of V4 on the momentum is a reflection of the fact that annihilation occurs, 
on an atomic scale, nearly at a point in coordinate space. At the level of precision required 
here we will also need to consider a modification that accounts for the interaction not being 
exactly pointlike [[LI]], 



Vr(k,l)= \ 7m / (k 2 + l 2 ), (8) 

where X = (19tt 2 - 132)/(tt 2 - 9). 

In first-order perturbation theory, taking the expectation value of V4 (which corresponds 
to multiplying by the square of the wave function at the origin m 3 a 3 /87r) and using T = 
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-2lm(E) reproduces Eqn. 2. In addition to this amplitude other operators accounting for 



relativistic effects are present ||10|| , and lead to the following ultraviolet divergent expression 
in second-order Rayleigh-Schrodinger perturbation theory, which has the first order effect of 
l/ 4 der included, 



ot ,ct,o \ r 8A 1, A 11 X ,2A 37T,-, 

nrqed = r 1 + -e x + (-) 2 e 2 ) + « 2 r [--— - -In- - - - — — - — . (9) 
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Here we have renormalized the imaginary part of V4 with a power series in a and introduced 
an ultraviolet cutoff mA on the momentum space integrations. 

The constants e\ and e 2 are determined by requiring that the amplitude for free particle 
scattering at threshold in NRQED be equal to that determined in a complete QED calcu- 
lation. The one-loop QED calculation at threshold has an amplitude corresponding to the 
decay rate 

ri = -r [- + A(A)] (10) 

7T A 

where A(X) = —10.28660 + 15.39A. Even though the limit A — > is taken at the end of the 
calculation, we keep terms of order A in the one-loop calculation because some terms enter 
the two-loop calculations with a factor 1/A. The one-loop matching calculation then allows 
us to determine 

8A XA im 7rrX XnX 

ei = T + — + AW ~l2~ ■ 

If we further define the two-loop QED decay rate as 

rH ? fri| ^! + ^ + ^ + S! l (12) 

the two-loop matching calculation gives 

ti\ n 2 X IItt 2 27r 2 ln2 „ 

e 2 = — InA + + B 2 . 13 
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The reason for defining T 2 in terms of A(X) rather than the physical limit A(0) is a practical 
one, associated with subtraction schemes required to deal with the most infrared divergent 
two-loop diagrams. It also has the advantage of leading to an exact cancellation with the 
factor A(X) in e\ in the matching calculation. 

With this determination of e\ and e 2 the NRQED decay rate becomes ultraviolet and 
infrared finite, and is given by 



a .. . .a. 9r 7r 2 , o/ll 21n2. 9 X 
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and numerically the constant we wish to determine is given by B = B 2 +tt 2 (0. 9129+1. 3302) = 
B 2 + 22.14. We note that the constant 0.9129 differs from Ref. fll |, where it is given as 1.16, 



the numerical value of 13/8-2/3 ln(2). An additional contribution of -1/4 to this number has 
recently been found [12[], removing the discrepancy. In addition Hill and Lepage [|l3j have 



recently recalculated a number of QED effects in a new nonperturbative implementation 
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of NRQED, and obtain 0.9125(5), so all NRQED calculations are now in agreement. As 
an additional check of the method, we verified that our implementation of NRQED, when 
applied to the one-photon annihilation contribution to ground state hyperfine splitting, 



reproduces the known answer |T0 



To finish the calculation the free two-loop QED calculation must be carried out and 
B2 extracted. While the two-loop calculation is sufficiently involved that we must defer 
its detailed description to a longer work [13], it is useful to refer to the diagrams that 



enter the one-loop calculation. The six diagrams of Fig. 1 we refer to as the outer vertex 
(OV), inner vertex (IV), double vertex (DV), self-energy (SE), ladder (L), and annihilation 
(A) contributions. After the ultraviolet divergences are removed by renormalization, the 
individual values of the diagrams are presented in Table I. 

The 81 QED diagrams that contribute to the decay rate at two- loop order break into 
11 classes that we label a-k. Class a consists of irreducible two-loop vertex corrections 
which generalize the OV diagrams. While free of binding singularities, their evaluation is 
complicated by the need to carry out two-loop renormalization. Similar comments apply to 
Class b, the two-loop generalization of the IV diagram, and Class c, the generalization of 
the SE diagram. We regulate ultraviolet divergences by using dimensional regularization, 
working in n = 4 — 2e dimensions, and the finite photon mass regulates infrared divergences. 
Renormalization constants in this scheme have not to our knowledge been presented in the 
literature: details of their calculation will be given elsewhere [ 15] . 

Class d consist of diagrams with reducible two-loop corrections, in which two separated 
ultraviolet divergent one-loop corrections are present. The next set of diagrams, which have 
no ultraviolet or infrared singularities, are those of class e, which generalize the DV diagrams. 
The most difficult to evaluate diagrams were in the f class, which generalize the L diagram. 
Most of these contributions diverge as 1/A, and the most singular as 1/A 2 . Canceling In A/ A 
divergences characteristic of Feynman gauge were present that were quite difficult to handle 
numerically. 

Class g consists of 9 diagrams in which the DV diagram has ultraviolet divergent one-loop 
radiative corrections in all possible vertices and propagators. Class h consists of radiative 
corrections to the A diagram, and have previously been calculated in Ref. [[HJ. Because that 
calculation used a Bethe-Salpeter formalism a lna was present that has to be replaced with 
a InA in our present formalism: the additive constant however is unchanged. Class i consists 
of diagrams where a vacuum polarization loop has been inserted in all possible places in the 
one-loop calculation. These have also been previously treated |T7| and |TB[ , as have the last 



two- loop effect we include, the square of the one- loop amplitude, which we call class j [19 



B. As with class h, our present formalism leads to terms that depend on A, but the additive 
constant is again unchanged. 

Finally, class k, which involves two of the three photons emitted in the decay undergoing 
light-by- light scattering have not been calculated: because of the small numerical contri- 
butions of these diagrams in parapositronium || we consider it highly unlikely that the 
omission of these diagrams will affect our conclusions. 

The results of the calculation are summarized in Table II. 
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TABLES 



TART F T 


1 < n-n rw m o i r7nri 


one-loop contributions to the orthopositronium 


LitJCdj Idle. 


contribution 


f(f) r o 


lnA(2)r„ 






Lov 





-4 


-1.028861425 


1.8756(1) 


Tiv 





-2 


-1.839322925 


4.7124 


Tdv 








-3.567629(21) 


7.5499(2) 


TSE 





4 


4.784983909 


-11.0445(1) 


r L 


2 


2 


-7.821768(32) 


12.296(4) 


r A 








-0.8140573(3) 


0.0 



total 2 -10.286606(10) 15.389(4) 



TABLE II. Contributions to the orthopositronium decay rate by class. 



diagram 


a 2 t 


a 2 -p 


lnA^ 

TT 


(f) 2 ln 2 AL 


(f) 2 lnAL 


(f) 2 r 


a 








-Lov 


-2 





-5.618 


b 








-Liv 


-1 





-0.705 


c 








-Lse 


2 





0.058 


d 

















2.421 


e 

















9.259(9) 


f 


2 In 2 


A 


LsE+OV+IV+DV 


1 


-2vr 2 /3 


-20.50(26) 


g 








— Ldv 








-1.372 


h 














TT 2 


9.007 


i 

















0.965 


j 


1 


A 











28.860(2) 


total 


21n 2+1 


2A 








7T 2 /3 


22.38(26) 
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We see from Table II that B 2 = 22.38(26), so the complete result for B is 44.52(26). 
While this is indeed a relatively large contribution, it leaves the theoretical prediction, 



T = 7.039 934(10)/zs~ 1 (15) 

well below the Ann Arbor results, by 8 and 5 standard deviations for the gas and vacuum 
measurements respectively, though consistent with the Tokyo measurement. 

No conclusions can be drawn until the experimental situation is clarified. If the Ann 
Arbor results are confirmed, while of course it is tempting to consider explaining this effect 
through exotic interactions ||20|| , it is worth noting that there is also at present a signifi- 
cant discrepancy between theory and experiment in the ground state hyperfine splitting of 
positronium. While there were disagreements between various calculations for some time, 
recently complete agreement ||22|| , [23[j, [fX3f] on a value of 203 392.05 MHz has been 



found. This value is 4 experimental standard deviations above the Yale measurement |24 
and 2.8 above the Brandeis measurement ||25|| . If the present disagreement of theory and 
experiment in positronium persists, any explanation in terms of new physics would be most 
compelling if both discrepancies were explained by the same mechanism. 

The work of JS was partially supported by NSF grant PHY-9870017, and that of GA by 
NSF grants PHY-9711991 and PHY-9722074. GA acknowledges the hospitality of the Aspen 
Center for Theoretical Physics and the Physics department of UCLA. We thank P. Labelle, 
G.P. Lepage and R. Hill for very useful conversations, and the constant encouragement of 
the Michigan experimental group is also gratefully acknowledged. 



6 



REFERENCES 



C.I. Westbrook, D.W. Gidley, R.S. Conti, and A. Rich, Phys. Rev. A 40, 5489 (1989). 

J.S. Nico, D.W. Gidley, A. Rich, and RW. Zitzewitz, Phys. Rev. Lett 65, 1344 (1990). 

S. Asai, S. Orito, and N. Shinohara, Phys. Lett. B 357, 475 (1995). 

W.E. Caswell, G.P Lepage, and J. Sapirstein, Phys. Rev. Lett. 38, 488 (1977). 

G.S. Adkins, Phys. Rev. Lett. 76, 4903 (1996). 

W.E. Caswell and G.P. Lepage, Phys. A 20, 36 (1979). 

S.G. Karshenboim, Zh. Eksp. Teor. Fiz. 103, 1105 (1993) [JETP 76, 541 (1993)]. 

A. Czarnecki, K. Melnikov, and A. Yelkhovsky, Phys. Rev. Lett. 83, 1135 (1999) and 

preprint |hep-ph/99104~88l . 



W.E. Caswell and G.P. Lepage, Phys. Lett. 167B, 437 (1986). 
A.H. Hoang, P. Labelle, and S.M. Zebarjad, |hep-ph/9909495| . 
P. Labelle, G.P. Lepage, and U. Magnea, Phys. Rev. Lett. 72, 2006 (1994). 
P. Labelle, private communication. 

R. Hill and G.P. Lepage, Cornell preprint (2000) and private communication. 
G.S. Adkins, R.N. Fell, and J. Sapirstein, in preparation. 
G.S. Adkins, R.N. Fell, and J. Sapirstein, in preparation. 
G.S. Adkins and M. Lymberopoulos, Phys. Rev. A 51, 2908 (1995). 
A.P Burichenko and D. Yu. Ivanov, Phys. At. Nucl. 58, 832 (1995) [Yad. Fiz. 58, 898 
(1995)]. 

G.S. Adkins and Y. Shiferaw, Phys. Rev. A 52, 2442 (1995). 
A.P. Burichenko, Phys. At. Nucl. 56, 640 (1993) [Yad. Fiz. 56, 123 (1993)]. 
A. Czarnecki and S.G. Karshenboim, [hep-ph / 99 1 14 1 0| . 
K. Pachucki, Phys. Rev. A 56, 297 (1997). 

A. Czarnecki, K. Melnikov, and A. Yelkhovsky, Phys. Rev. Lett. 82, 311 (1999); Phys. 
Rev. A 59, 4316 (1999). 

G. Adkins and J. Sapirstein, erratum to Phys. Rev. A 58, 3552 (1998). 
MW. Ritter, P.O. Egan, VW. Hughes, and K.A. Woodle, Phys. Rev. A 30, 1331 (1984). 
A.P. Mills, Jr., Phys. Rev. A 27, 262 (1983). 

Figure captions: 

Fig. 1 Diagrams contributing to the one-loop decay rate correction. 
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FIG. 1 



